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Abstract 

oo i Risken's matrix continued fraction method is used to solve the Fokker-Planck equation to calcu- 

late particle current in an inertial symmetric (sinusoidal) periodic potential under the action of a 

(—1 '. 

Q . constant force. The particle moves in a medium with friction coefficient also varying periodically in 

(P . 

space as the potential but with a finite phase difference (j){j^ mr, n = 0, 1, 2, ...). The algebraic sum 

of current with applied forces ±1^1 gives the ratchet current in the adiabatic limit. Even though, 

■^ ■ the effect of frictional inhomogeneity is weak, the ratchet current shows very rich qualitative char- 

acteristics. The effects of variation of F, the temperature T, and the average friction coefficient 70 

__i ^ on the performance characteristics of the ratchet are presented. 
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I. INTRODUCTION 

In equilibrium, it is not possible to get a steady current in a periodic potential without 
the application of external forces. This is in conformity with the second law of thermody- 
namics. However, in recent times it has become a popular theme to obtain such currents in 
nonequilibrium situations. The currents so obtained in asymmetric but periodic potentials 
with no imposed directional preferences in the presence of thermal noise are called thermal 
ratchet currents jj, l2]- The usual methods one generally employs to realise such currents are: 
Rock the system periodically alternately with constant forces |F| and —\F\ for equal short 
intervals 3|, or flash the potential eitherperiodically or randomly with two unequal (for 
example, zero and nonzero) amplitudes [J]. However, the method of obtaining current in 
symmetric periodic potentials subject to asymmetric periodic forcings with zero mean per 
period has also been established 5|. In most of the cases the particle motion is considered 

)lecular motor) transport 



61] . There is yet another 



overdamped which is more suitable to model particle (such as a mo 

in biological systems where large friction plays important role [sl, 14, 

model wherein the friction is considered nonuniform in space and in particular periodically 

varying similar to the potential function. The ratchet current could be obtained, in this 

inhomogeneous case, even in symmetric periodic potentials driven periodically and sym- 

metncally ■„ t,me0. So far such currents are obtained „r overdarnpedH mhonrogeneous 

systems and the possibility of obtaining such currents in underdamped systems [91] had only 



been indicated earlier. This model is very appropriate to describe Josephson junctions [10 ] 
(where the frictional inhomogeneity bears an exact analogy to the coupling between con- 
comitant Cooper pair and the quasiparticle tunnelings), and motion of adatoms on crystal 



surfaces llj . In this part (I) of the work we present detailed characteristics of the ratchet 
current in inhomogeneous (nonuniform friction) inertial systems when driven adiabatically. 
And in the next part (II) we shall present results for similar systems when driven at finite 
frequency. 

As will be described in the next section, the steady state mean velocity v of a particle 
under the action of a constant force F is calculated using the matrix continued fraction 
method (MCFM). The algebraic sum S-u (= 'y(]i^]) + v{—\F\)) of these mean velocities v 
with applied forces \F\ and —\F\ is termed here as the ratchet current in the adiabatic limit. 
Since the phase difference (p{y^ 0, nn) between the sinusoidal potential V{x) and the friction 



coefficient 7(x) is the only parameter in the model that breaks the directional symmetry of 
the problem, the ratchet current is expected to be small. However, Sw turns out to be an 
interesting complicated function of F, the temperature T, and the mean friction coefficient 
7o. This is not surprising and can be understood quite simply on physical grounds as 
discussed elsewhere [l2f|. 
In the same work 



12| | the variation of the ratchet current "Ev as a function of applied 
force \F\ for various values of the average friction coefficient 70 at temperature T = 0.4 has 
been presented. (The same choice of dimensionless parameters, T, etc, are elaborated in 
Sec. II below.) It shows that at the smallest explored 7o(= 0.035), Sv peaks, with value 
^ —0.08, close to F = 0.09 ^ Fci= 1); the critical field. The current is thus thermally 
assisted. The variation of Sv is very sharp: Appreciable current could be seen only in the 
range 0.05 < |F| < 0.12. The current strength reduces almost to zero for F = 0.12. Beyond 
this F no meaningful v could be obtained from MCFM. As 70 is increased the range of F, 
where appreciable T^v results, widens monotonically. However, the peak current becomes 
maximum for 70 ~ 0.4, with maximum current c^ —0.28 at |F| c^ 0.75, and then gradually 
decreases as 70 is increased further. By a close observation of the reported figure (3, ref. 



12l |) one can conclude that Sw at T = 0.4 will show nonmonotonic variation with 70 at 
carefully chosen F values, which as we present in the following is, indeed, the case. 

For a given applied force F, and the temperature T the current S^ initially increases, 
attains a peak, and then gradually decreases as 70 is increased. The result is reminiscent 



of the well-known result [13|] obtained by Kramers, however, in an entirely different context 
of particle flux across a potential barrier. In all situations T,v increases with \F\, attains a 
peak, and gradually decreases. With appropriately chosen parameters 70 and T, T,v shows 
current reversal as well and consequently as |F| is increased further St; changes its direction 
and increases in magnitude eventually again decreases in magnitude after attaining a second 
peak. Interestingly, Ev shows very similar behaviour also as a function of temperature T in 
some parts of the parameter space (70, F). 

As the value (and direction) of S-u changes one naturally asks about how it is correlated 
with the changes in energies of the system. It turns out that the average kinetic energy of 
the particle increases monotonically with \F\, whereas the average potential energy shows a 
nonmonotonic behaviour. Again, in this case too we consider the mean of the average kinetic 
(as also of the potential) energies for F = ±|F|. The average kinetic energy overwhelms 



the potential energy as F is increased forcing the average total energy also to increase 
monotonically. However, the contribution of the average potential energy of the particle is 
reflected in the derivative of the average total energy with respect to F. These derivatives 
show distinctive behaviour where S^ peaks as a function of F. 

In order to be useful for practical purposes the system must have good performance 
characteristics. We calculate the efficiency of energy transduction of this inhomogeneous 
inertial system driven adiabatically and symmetrically and moving in a symmetric periodic 
potential. The system yields small current and hence the work done against an applied load 
is also small if at all it is possible to obtain. Naturally, the thermodynamic efficiency of the 
system is also small and it is feasible to obtain only in a very small region of parameter space. 
However, even without an applied load the system continuously performs work against the 
frictional force and in order to take it into account a quantity called Stokes efficiency has 



been defined |14j. Il5l . Il6| . We calculate the Stokes efficiency of our system. This efficiency 
also has a complex behaviour and, in particular, it is nonmonotonic. It is also possible, 
interestingly, to observe increasing Stokes efficiency with temperature contrary to what is 
observed in usual macroscopic heat engines. This kind of behaviour has been shown to occur 



in overdamped inhomogeneous but asymmetrically driven systems 171 ]. 

Admittedly, the quantitative effect of the inhomogeneity is not very large. Nevertheless, 
our investigation shows that many important results can be obtained, some even counterin- 
tuitive, albeit at a small but nonnegligible scale, using system inhomogeneity. The ratchet 
current, the efficiency of performance, and so on, can be improved by orders of magnitude, 
if zero mean asymmetric periodic drive is employed instead of a symmetric drive. In this 
case, however, inhomogeneity does not play a decisive role. Similar results can be obtained 
without system inhomogeneity as the asymmetric drive plays the major role. 

The above mentioned results will be elaborated upon in section HI. The paper will be 
concluded with a discussion in section IV. 

II. THE MATRIX CONTINUED FRACTION METHOD FOR INHOMOGE- 
NEOUS INERTIAL SYSTEMS 

The motion of a particle of mass m, moving in a periodic potential V{x) = —VQsin{kx) 
in a medium with friction coefficient ^{x) = 7o(l — \sin{kx + (;/>)) (0 < A < 1) and subjected 



to a constant uniform force F, is described by the Langevin equation JTI. 1 18l|. 

.^--7W$-^ + ^+V^WT^W. (2.1) 



dt'^ ' ^ ^ dt dx 



Here T is the temperature in units of the Boltzmann constant A;^. The Gaussian distributed 
fluctuating forces ^(t) satisfy the statistics: < ^{t) >= 0, and < C,(t)^{t ) >= 26{t — t). For 
covenience, we write down Eq.(2.1) in dimensionless units by setting m = 1, Vq = 1, k = 1 



so that, for example, the scaled 70 is expressed in terms of the old y/mVok, etc. The reduced 
Langevin equation, with reduced variables denoted again by the same symbols, is written as 

d X dx I 

-^ = -7(3^)^ + cosx + F+ ^^{x)Ti{t), (2.2) 

where 7(x) = 70(1 — \sin[x + (j))). Thus the periodicity of the potential V{x) and also 
the friction coefficient 7 is 27r. The barrier height between any two consecutive wells of 
V{x) persists for all F < 1 and it just disappears at the critical field value F^. = 1. The 
noise variable, with the same symbol ^, satisfies exactly similar statistics as earlier. The 
Fokker-Planck equation corresponding to Eq.(2.2) is given as 

dW{x,v,t) 
di 

where the Fokker-Planck operator 



CFpW{x,v,t), (2.3) 



d d d 

Cfp = -v^r + 7o(l - \sin{x + (j)))Tr'^ ~ {cosx + F)—- + 
ox ov ov 

7oT(l-As^n(x + 0)) — , (2.4) 

and W{x, v, t) is the probablity distribution of finding the particle at position x with velocity 
f = ^ at time t. 

It is hard to solve the equation in the general case. However, our aim is to solve Eqs. (2.3- 

n 

4) in the steady state case. We apply MCFM developed by Risken, et. al. [19l], and used by 



others to suit various circumstances [20|. The basic idea behind the method is to expand the 
distribution function W{x,v,t) in series in terms of products of (Hermite) functions i^niv), 
which are explicit functions of velocity v, and Cn{x,t) (as coefficients) which are explicit 
functions of position x and time t: 

-1 - 2 °° 
W{x, V, t) = (27rr)-e^ ^ C„(x, t)^n{v), (2.5) 

n=0 



with 



ijn = —T=^, (2.6) 



and 

i)o = {27v)^T^e^. (2.7) 

It turns out that ^„(f ) (n = 0, 1, 2, ...) are the eigenf unctions of the operator — 7(x)6^6: 

- -f{x)b^btpn{v) = -'j{x)niJn{v), (2.8) 

where b and b^ are the bosonic-hke annihilation and creation operators satisfying 

66t _ 5t^ = i^ (2.9) 

and defined as 

(Note: b^ipn{v) = y/n + l^pn+i{v) and bip^iv) = y/nijn~i{v) ■) 

Substituting (2.5) in the Fokker-Planck equation (2.3) and multiplying from left by e^T 
on both sides of Eq. (2.3) and on simplifying, we get, 

f; ^^'"j^^'^V n. = Cfp f; C^(X, t)i^m{v), (2.12) 

where 

CFP = e^CFpe^. (2.13) 

The operator Cpp is given by 

£FP = -^{xpb-bD-b^D, (2.14) 

where the operators D and D are defined as, 

D = Vf-^, (2.15) 

fj T/' p 

with the prime on the potential V denoting differentiation with respect to x. 



Multiplying both sides of Eq.(2.12) from left by ipn and integrating over the velocity 
variable and using the orthonormality properties of the functions ^„, we get the equation 
on Cn{x,t): 

^^J^^A + n7(x)C„(x,t) + V^^TlDCn+i{x,t) + v^DC„_i(x,t) = 0. (2.17) 

In the stationary (steady) state case, we thus have, 



n-i{x)Cn{x) + y/^r+lDCn+i{x) + V^DCn-i{x) = (2.18) 

giving a recursion relation on C„(x) (n = 0,1,2,...). The first (n = 0) of the recursion 
relations reveals that Ci is a constant. 

In this steady state case, the normalization condition is given by 

/»2-7r /»oo /»2-7r 

/ / W{x,v)dxdv= Codx = l, (2.19) 

Jo J~oo Jo 

giving a condition on the coefficient Co{x). The mean velocity v{F), the average potential 
energy Epot, and the average kinetic energy E^in, are obtained, respectively, as, 

/'27r /»oo 

v= / vW{x,v)dxdv = 2TTVfCi, (2.20) 

Jo J-oo 

/»7r 

Epot= / V{x)Co{x)dx, (2.21) 



and 

^2 1 ^ .2. 



^fc^n =< y >= 2 + 7f / '^2(a:)(/a;. (2.22) 

Thus, all the physical quantities of our interest can be calculated provided Co, Ci, and C2 
are known from Eq.(2.18). 

Writing Cn{x) = EgQe^«^. V{x) = Eg'^"^e^«^. and 7(0;) = ^^ ^5e*'?(^+'^) , where i = 
^/—l and g's as integer indices the recurrence relation (2.18) becomes. 



p p p 

-^ E F^^,pC:_, + ^ E ^(^ - P)^'^"'C'L,. (2.23) 



(7)^'^ 


' = !'- 


PQi{Q-p)<t> 




(D)«'P 


= ipSq,p, 




(1)^ 


hP — X 


y'y,p 


= i{q- 


-p)V^-''., 



Defining the matrices 7, D, 1, and V with components 

(2.24) 
(2.25) 
(2.26) 
(2.27) 

and the column matrices Cn, n = 0, 1, 2, .. with components 

(C„)^ = CI (2.28) 

the recurrence relation can be written as, 

n7Cn + v/^7TTv^DCn+i + v^(v^D + ^[V' - Fl])Cn_i = 0. (2.29) 



Introducing the matrices Mn, (n = 0, 1, 2, ..) such that 

v^DCn = M„Cn_i (2.30) 

the recursion relation (2.29) for n = 2 becomes 



M2Ci(= V2DC2) = -D(7 + ^M3)-^(v^D + ^[V' - Fl])Ci. 

2 vT 

Hence, for nonzero Ci, M2 can be given in terms of M3 as 



(2.31) 



M2 = D(7 + ^M3)-^(-^[Fl - V'] - v^D). (2.32) 



and similarly, one can write down for M3, M4, and so on. And, in general, 

M„ = D(7 + ^M„+i)-^(^[Fl - V'] - v^D), (2.33) 

n VT 

for 2<?7,<A^ — 1, if the series of recursion relations is terminated at n = A^ (Cn = for 
n> N + 1). Thus, the last {n = N) relation becomes, 

Mn = T>-f~\^[Fl - V'] - VtD) , (2.34) 



for N > 2. Hence, Mn-i, Mn-2, ••• M2 can be given entirely in terms of 7, V , and Fl. 
Using v^DC2 = M2C1, the recursion relation (n = 1) 



7Ci + V2VTBC2 + {VTB + ^[V - Fl])Co = 
yields 

Co = HCi, (2.35) 

where 

U={^[Fl-Y']-Vf-D)-\^ + VfM2). (2.36) 

vT 



Since Ci is a constant, 

C7« = ^(H)«'P(Ci)P = (H)«'°Ci. (2.37) 

p 

From the normahzation condition (2.19), we obtain 27r(Co)° = 1, and hence from Eq.(2.36), 



1 

27r(H) 



Ci = TTZT^^o- (2-38) 



Therefore, once H is calculated Ci is known and hence all the components Cq = (Cq)^ are 
evaluated. 

III. NUMERICAL RESULTS 

In order to calculate the physical quantities of interest, namely, the mean velocity v, the 
average kinetic and potential energies for given values of F, 70, and temperature T one needs 
to evaluate the matrix H. For our special case, the potential V{x) and 7(x) each has only 
two Fourier components. Therefore, the dimension of H depends on the number Q of Fourier 
components to be kept for the coefficients C„(x). Also, the degree of complication increases 
with the number A^ of recursion relations one requires to be considered. The acceptable 
solution will require the matrix H to remain essentially unchanged as Q and A^ are changed. 
All these numbers depend on the values of F, 70, and the temperature T. Typically, in 
the worst case, it is sufficient if Q < 30 and A^ < 2000 is taken for small 70 <^ 1 in the 
appropriate parameter space range of {F,T). Outside this range no acceptable solution is 
found. It turned out that the paramater space {'~fo,F,T) range where solutions could be 
extracted was just the region that yielded finite ratchet current S-u. In all our calculations 
for which the results are described in the following we have set (p = 0.35 and A = 0.9. 







\ ~ \ F . 0.66 














\ \ F=1.0 


^^^^-^ 










-0.05 


] ' 


'^"^ 




^, 







-0.1 


'' / 


,.,- 


-■'■" 




" 


5 


-o.ts 


1 ' / 

\/ 

X 


/ 


-0.04 


F = 0.1 


- 




-0.2 


/\ ./ 


& 


-0.06 




- 




-0.25 










- 




0.044 








, 




To 





0.2 0.4 0.6 o.a 



1 1.2 1.4 1.6 1.! 

To 



FIG. 1: Shows the ratchet current Sw as a function of 70 for two values of F . The inset shows the 
variation of current vs 70 for F = 0.1. 



-0.02 y T = 0.6 
T = 2.0 
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FIG. 2: Ratchet current versus F is shown for 70 = 0.1 and two values of temperature. The inset 
shows variation of TiV for smaller values of -F at T =0.6. The dotted line for zero current is shown 
just to guide the eye. 

A. The ratchet current 

Figure 1 shows how Hv changes with 70 at T = 0.4. Contrary to usual perception current 
increases as the system becomes more damped. This trend is seen for lower 70 values. 
However, in the higher range of 70, as expected, the current diminishes with damping. 
v(|F|) and v(— |F|) are both monotonically decreasing in magnitude with 79. It is only their 
algebraic sum Sv that shows the nonmonotonic behaviour. It is purely an effect brought 
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rt-|F|) ■ 




0.2 0.4 



FIG. 3: Shows the variation oi ii(\F\) and iJ,{—\F\) with F for 70 = 0.5. The Uttle vertical marking 
(I) indicates the position of current maximum. 




FIG. 4: The variation of ratchet current is shown for 70 = 0.035 for different values of F. T,v 
decreases slowly with T(> 4). The inset shows two current reversals for 70 = 0.1, F = 0.1. The 
zero T,v line is shown to guide the eye. 

about by the symmetry breaking due to nonzero phase difference (p and any analogy with 
potential barrier crossing phenomenon at low damping described by Kramers is superfluous. 
In Fig. 2, the variation of Sv as a function of |F| for 70 = 0.1 is shown for two values 
of temperature T = 0.6 and 2.0. It is to be noted that T = 2.0 corresponds to an energy 
exactly equal to the potential barrier height at F = 0. One can see that in the scale of 
the figure the current Sv remains negative for T = 0.6, whereas at T = 2.0 the current is 



11 



positive at low |F| but it changes sign and remains negative at higher \F\ values for the 
same phase lag (p. This current reversal, however, is neither confined to higher T nor to 
higher 70 values. The inset of the figure shows that at smaller |F| region S-u is positive and 
then it becomes negative as \F\ is increased. Of course, in this case the positive current is 
an order of magnitude smaller and hence is not visible in the larger figure. A similar result 
is obtained for 70 = 0.035 where again the negative current at larger |F| dominate but 
positve currents at low \F\ region are also appreciable. Here again w(|-F|) and ^;( — |F|) both 
increase monotonically with \F\ only their relative values change giving rise to nonmonotonic 
behaviour of S-u, including its direction reversal. At higher values of \F\ the current should 
saturate with a small value since the mobility /i = |^| for both ±F tend to 7q"^. The 
variation of /i as a function of |F| is shown in Fig. 3. It shows that at the position of the 
current maximum /u(|-F|) = /i(— |-F|) as it should be. 

The variation of St; as a function of temperature T at various values of applied force 
\F\, as shown in Fig.4, reveals that current can also increase with temperature. Moreover, 
one can achieve current reversals as well. In Fig.4, plots are given for 70 = 0.035 and 
\F\ = 0.025, and 0.075. |F| = 0.075 yields large currents and shows one current reversal. 
For larger |F| the current remains negative for all T. This is not shown here to preserve 
clarity of the figure. However, in order to highlight the existence of two current reversals 
we plot the graph for 70 = 0.1, and \F\ = 0.1 in the inset of the figure. For low values 
of \F\, for example \F\ = 0.025, the current always remains positive, though small, which 
is consistent with the observation made with reference to Fig. 2. It is noteworthy that the 



magnitude of current obtained here is quite large compared to what has been reported 171] 
for overdamped systems. 

B. The energetics 

The average potential and kinetic energies can be easily calculated. The expressions are 
given in the Appendix. Fig.5 shows the variation of the average kinetic (Ekin), potential 
{Epot)i and total (-Etot) energies at T = 0.4 for 70 = 0.4 as a function of \F\. Monotonically 
increasing nature of Ekm-, and Etot and nonmonotonic behaviour of Epot are seen for all 
(explored) 70 and T values. Interestingly, the slope of Etot shows a shoulder peak which 
after reaching a minimum rises again. Fig. 6. For a given temperature the shoulder peak 
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FIG. 5: Shows the variation of the average kinetic, potential and total energies at T 
7o = 0.4 as a function of F. 
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FIG. 6: In fig. [a] the slopes of the average total energies are shown for various values of 70 < 1.0. 
The graphs for 70 = 0.2, 0.4, 0.5 and 0.75 are scaled by factors of 10, 20, 25 and 30, respectively, 
to make them visible in the figure with the graph for 70 =0.1. Fig. [6] shows plots for 70 > 1.0. The 
slopes for 70 = 2.0, 4.0 and 5.0 are blown up by 4, 5 and 5 times, respectively, to scale with the 
curve for 70 =1.0. The little vertical markings on the curves show positions of the current maxima. 

occurs for small 70 and large 70 values. However, there is an intermediate range of 70 values 
where the slope does not show this peaking behaviour; the total energy does not acquire 
negative curvature in this range of 70. The slope gives valuable information about the 
variation of the ratchet current St; as shown in Fig. 6. For small 70 < 1 values the slope 
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Region with negative curvature 
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FIG. 7: T — 7o plane showing two disjointed regions with possibihty of negative curvatures of the 
total energy Etot{\F\) separated by a region without the negative curvature of Etot{\F\)- 
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FIG. 8: Stokes efficiency versus F for different values of 70. The inset shows ry^tofces and Su plotted 
for 7o = 0.75. 

peaks at an intermediate value of |F|, and this is precisely where the current maximum 
occurs. For 70 > 1 the positions of the current maximum and that of the slope of the total 
energy are clearly well separated. In fact, the current maximum occurs, for these large 70 
values, even beyond the position of the minimum (after the peak) of the slope of the total 
energy. At this point it is worth mentioning that for not too large temperatures for 70 < 1 
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FIG. 9: Shows the variation of Stokes efficiency with T for various values of F for 70 = 0.4. The 
Uttle vertical markings on the efficiency curves indicate the positions of the current maxima. The 
inset shows efficiency and current for F = 0.7 

the current peak occurs at |F| < 1 whereas for all 70 > 1 the current peaks invariably 
occur at |F| > 1. Recall that |F| = 1 is the critical field at which the potential barrier to 
motion just disappears. A careful observation of the results at various temperatures T, thus, 
allows one to construct a diagram in the (T — 70) plane, Fig. 7, with two disjointed regions 
with possibility of negative curvatures of Etot{\Fo\) separated by a region with only positive 
curvatures. In the " negative curvature" region with 70 < 1 the slope of Epot peaks close to 
where Sv becomes maximum. Whereas in the other negative curvature region with 70 > 1 
current maximum occurs far away from the position (Fq) where the slope of Etoti\Fo\) peaks. 
It shows that one can clearly separate low damping (70 < 1) region from the high damping 
(70 > 1) case with each showing qualitatively different characteristic behaviour of particle 
motion. 

C. The efficiency of performance 

The currents obtained by application of symmetric drives ±|-F| are small and not sustain- 
able against applied loads of perceptible magnitude. Therefore it is not feasible to calculate 
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work done by the system against the load and hence it is difficult to study its behaviour 
as the other parameters vary. However, the particle keeps moving, essentially with average 
velocity Sv, against the frictional force and hence performs some work, though not useful 
practically. How efficiently the system does this work is obtained by comparing the mag- 
nitude of this work (= 7o(X]^)^) with the input energy supplied to it. A measure of this 



quantity called the Stokes efficiency is defined and given (in our case) asjl4l. llsLIld Il7l|. 



^' |F||(z;(|F|)-.;(-|F|))r ^'^ 

In Fig. 8, we plot r]s as a function of |F| for various 70 values at T = 0.4. r]s behaves in 
a very similar fashion as the corresponding current does. However, rjs and current do not 
peak at the same value of F as will be clear from the inset of the figure where r]s and Su 
are plotted together for 70 = 0.75. It is remarkable that with such a modest approach one 
can find regions in the parameter space (70, |-F|) where the efficiency of performance exceeds 
two per cent. It shows that exploiting system inhomogeneity to obtain ratchet current in 
inertial systems may be useful for practical purposes too. 

The variation of Stokes efficiency 775 with T is shown in Fig. 9 for \F\ = 0.66,0.70, and 
0.80 for which the efficiency and current are appreciable. It is noteworthy that at low 
temperatures the performance of the system improves with rise in temperature and only 
after attaining a peak it declines with further rise in temperature. In this case, however, the 
efficiency and current seem to peak at the same temperature as can be seen from the inset 
of the figure. For F = 0.66 we could not see the low temperature behaviour because we do 
not get sensible results in this region from MCFM. 

D. The asymmetrically driven ratchet 

It is clear that in general the efficiency of performance of the inhomogeneous inertial 
ratchets is small except in some special sections of the parameter space (70, |-F|,r). This 
efficiency can be dramatically enhanced if the ratchet is driven asymmetrically with an 
asymmetry parameter e(0 < e < 1). Here the driving force in the positive and in the 
negative directions are not same: The force F = jz^\F\ is applied in the positive direction 
for a fraction of time r = ^ whereas F = — |F| is applied for a fraction of time 1 — r = ^^. 
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F. 0.017 , T=0.035 - 
F = 0.05 , 1=0.1 
F = 0.27 . 7=0.4 - 




0.1 0.2 0.3 0.4 



FIG. 10: Stokes efficiency as a function of the asymmetry parameter e for different combinations 
of values of F and 70. 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 



FIG. 11: The thermodynamic efficiency rjffi versus load for two values of F with fixed e = 0.5, 70 
0.035 and T = 0.4. 

The ratchet current now will be given by 



T,v = v^ + v' 



(3.2) 



with 



f \F\), 

1 + e 



M-m 



(3.3) 
(3.4) 
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The Stokes efficiency of the thermal ratchet motion is similarly defined as in Eq. (3.1): 

Vs = |p||_,n+,x —T- 3.5 

In this case the ratchet current, in a typical situation, is quite large. Correspondingly, the 
Stokes efficiency also become large. Fig. 10 shows the variation of r]s as a function of the 
asymmetry parameter e for various values of (a) \F\ =0.017, 70 =0.035; (b) \F\ =0.05, 
7o =0.1 and (c) \F\ =0.27, 70 =0.4 at T = 0.4. Clearly, more than 10 per cent of ratchet 
operation efficiency becomes achievable. Since the currents are large one can also drive these 
ratchet currents against applied loads, L. Hence, the thermodynamic efficiency rjth, defined 
in this case as 15|, UM, llTj , 

Vth = n^||_, .i+eN —J 3.6 

III V 1 — e' I 

can be calculated. The variation of rjth as a function of L for 70 = 0.035 at T = 0.4 is 
shown in Fig. 11 for \F\ = 0.03, and 0.053. From the figure one can see that thermodynamic 
efficiency can be substantial. One can find an optimum load, L, and an optimum asymmetry 
parameter, e, of field drive where the ratchet performs most efficiently. It should, however, 
be noted that the effect of asymmetric drive is overwhelming and the system inhomogeneity 
plays only a secondary role in this case. One can achieve almost similar results without 



invoking system inhomogeneity as hs been seen in the overdamped case 17| 



IV. DISCUSSION AND CONCLUSION 

All the results of inertial particle motion, in a symmetric periodic potential and driven 
symmetrically ( and of course adiabatically) , such as ratchet current and its nonmonotonic 
behaviour (including its direction reversals) as a function of \F\, T, and 70 are possible here 
only because the friction was considered space-dependent. The periodic variation of the fric- 
tion in space is just a model and we have not provided any microscopic basis for it. However, 
it has been shown earlier that in order to describe the motion of an adsorbed particle on a 
crystalline surface (of same atoms as the adsorbed one) by a Fokker-Planck equation suit- 
able for one-dimensional particle motion, one necessarily have to consider space-dependent 
friction HI] . It turned out from their microscopic theory that the space dependence has 
to be periodic as the average potential created by the crystalline (surface) substrate for 
the adatom and in antiphase (0 = it) with it. In other words, the friction will be large 
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where the substrate potential is small. The possiblity of obtaining current due to diffusion 



inhomogeneity 2l|] and frictional inhomogeneity 18| has been argued earlier. The phase dif- 
ference 7^ nn{n = 0, 1,2, ...) is necessary to break the spatial symmetry in order to get 
ratchet current. We have found that the ratchet current depends on the value of (p. For 
instance, taking = 0.4 gives different and larger current from what we have reported here 
for (p = 0.35 at the same point in the parameter space of (|F|,r, 70). 

The increase of ratchet current and also the efficiency of the ratchet as a function of T at 
small temperature range is because the particle has a large potential to surmount at small 
\F\ and increasing temperature helps in achieving that. However at the larger temperature 
range obtaining a second current peak as shown in Fig. 4 cannot be explained as simply and 
it arises because of the competition between the temperature and friction effect as argued 



elsewhere |l 21]. 



The behaviour of the slope of the total energy vis-a-vis the position of ratchet current peak 
shows a clear demarcation occurring around 70 = 1, in the (T — 70) plane. The investigation 
of time scales of running and locked states of the particle motion may provide a clue to this 
remarkable behaviour. We are carrying out numerical Langevin dynamic simulation on this 
aspect of the problem and the results will be reported elsewhere. We conclude by reiterating 
that particle motion in inhomogeneous systems has the potential of a rich field of research. 
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APPENDIX 

Using Eqs. (2.36) and (2.37), the potential energy expression (2.21) is given by 

:h)«'0 



E. 






q 



The kinetic energy expression (2.22) can be written as 

Ekin = 1 + v^7r(C2)°. 
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From Eqs. (2.29) and (2.32) one can write down a general equation relating Cn and 
Cn-i(^ > 1) as 

'n 



where the matrices Ani(r?7, = 0, 1, 2, ...) are given by 

A„, = (7 + Vf{m + l)-^'DK^+^Y\^[Fl ^ V'] - VtD). 
If the recursion is terminated at m = A^ (An+i = 0) 



for N > 2. Co is thus calculated as 



C2 — ^=AiAoCo. 
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